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Abstract: We develop a pricing rule for life insurance under stochastic mortality in an 
incomplete market by assuming that the insurance company requires compensation for its 
risk in the form of a pre-specified instantaneous Sharpe ratio. Our valuation formula satis- 
fies a number of desirable properties, many of which it shares with the standard deviation 
premium principle. The major result of the paper is that the price per contract solves 
a linear partial differential equation as the number of contracts approaches infinity. One 
can interpret the limiting price as an expectation with respect to an equivalent martin- 
gale measure. Another important result is that if the hazard rate is stochastic, then the 
risk-adjusted premium is greater than the net premium, even as the number of contracts 
approaches infinity. We present a numerical example to illustrate our results, along with 
the corresponding algorithms. 
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1. Introduction 

We propose a pricing rule for life insurance when interest rates and mortality rates are 
stochastic by applying the method developed and expounded upon by Milevsky, Promislow, 
and Young (2005, 2007). In the case addressed in this paper, their method amounts to 
targeting a pre-specified Sharpe ratio for a portfolio of bonds that optimally hedges the 
life insurance, albeit only partially. 

Actuaries often assume that one can eliminate the uncertainty associated with mortal- 
ity by selling a large number of insurance contracts. This assumption is valid if the force of 



1 



mortality is deterministic. Indeed, if the insurer sells enough contracts, then the deviation 
of actual results from what is expected goes to zero, so the risk is diversifiable. However, 
because the insurer can only sell a finite number of insurance policies, it is impossible to 
eliminate the risk that experience will differ from what is expected. The risk associated 
with selling a finite number of insurance contracts is what we call the finite portfolio risk. 

On the other hand, if the force of mortality for a population is stochastic, then there 
is a systematic (that is, common) risk that cannot be eliminated by selling more policies. 
We call this risk the stochastic mortality risk, a special case of stochastic parameter risk. 
Even as the insurer sells an arbitrarily large number of contracts, the systematic stochastic 
mortality risk remains. 

We argue that mortality is uncertain and that this uncertainty is correlated across 
individuals in a population - mostly due to medical breakthroughs or environmental factors 
that affect the entire population. For example, if there is a positive probability that medical 
science will find a cure for cancer during the next thirty years, this will influence aggregate 
mortality patterns. Biffis (2005), Schrager (2006), Dahl (2004), as weU as Milevsky and 
Promislow (2001), use diffusion processes to model the force of mortality, as we do in 
this paper. One could model catastrophic events that affect mortality widely, such as 
epidemics, by allowing for random jumps in the force of mortality. Also see the related 
work of Cox and Lin (2004) and Cairns, Blake, and Dowd (2004). 

In related work, Blanchett-Scalliet, El Karoui, and Martellini (2005) value assets that 
mature at a random time by using the principle of no arbitrage by focusing on equivalent 
martingale measures; the resulting pricing rules are, therefore, linear. Dahl and M0ller 
(2006) take a similar approach in their work. However, for insurance markets, one cannot 
assert that no arbitrage holds, so we use a different method to value life insurance contracts 
and our resulting pricing rule is non-linear, except in the limit. 

We value life insurance by assuming that the insurance company is compensated for its 
risk via the so-called instantaneous Sharpe ratio of a suitably-defined portfolio. Specifically, 
we assume that the insurance company picks a target ratio of expected excess return to 
standard deviation, denoted by a, and then determines a price for a life insurance contract 
that yields the given a for the corresponding portfolio. One might call a the market price 
of mortality risk, but it appears in the pricing equation in a non-linear manner. However, 
as the number of life insurance policies increases to infinity, then this a is the market price 
of risk in the "traditional" sense of pricing in financial markets in that it acts to modify 
the drift of the hazard rate process. 

In this paper, we assume that the life insurance company does not sell annuities to 
hedge the stochastic mortality risk. In related work, Bayraktar and Young (2007a) allow 
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the insurer to hedge its risk partially by selling pure endowments to individuals whose 
stochastic mortality is correlated with that of the buyers of life insurance. 

We obtain a number of results from our methodology that one expects within the 
context of insurance. For example, we prove that if the hazard rate is deterministic, then 
as the number of contracts approaches infinity, the price of life insurance converges to 
the net premium under the physical probability; see Corollary 4.4. In other words, if the 
stochastic mortality risk is not present, then the price for a large number of life insurance 
policies reflects this and reduces to the usual expected value pricing rule in the limit. 

An important theorem of this paper is that as the number of contracts approaches 
infinity, the limiting price per risk solves a linear partial differential equation and can 
be represented as an expectation with respect to an equivalent martingale measure as 
in Blanchett-Scalliet, El Karoui, and Martellini (2005); see Corollary 4.6. Therefore, we 
obtain their results as a limiting case of ours. Moreover, if the hazard rate is stochastic, 
then the value of the life insurance contract is greater than the net premium, even as the 
number of contracts approaches infinity; see Corollary 4.5. 

The remainder of this paper is organized as follows. In Section 2, we present our 
financial market, describe how to use the instantaneous Sharpe ratio to price life insurance 
payable at the moment of death, and derive the resulting partial differential equation 
(pde) that the price A = A^^^ solves. We also present the pde for the price yl(^) of n 
conditionally independent and identically distributed life insurance risks. In Section 3, we 
study properties of ^^"); our valuation operator is subadditive and satisfies a number of 
other appealing properties. In Section 4, we find the limiting value of ■^A^'^^ and show that 
it solves a linear pde. We also decompose the risk charge for a portfolio of life insurance 
policies into a systematic component (due to uncertain aggregate mortality) and a non- 
systematic component (due to insuring a finite number of policies); see equation (4.18). 
In Section 5, we present a numerical example that illustrates our results, along with the 
corresponding algorithms that we use in the computation. Section 6 concludes the paper. 

2. Instantaneous Sharpe Ratio 

In this section, we describe the term life insurance policy and present the financial 
market in which the issuer of this contract invests. We obtain the hedging strategy for 
the issuer of the life insurance. We describe how to use the instantaneous Sharpe ratio to 
price the insurance policy, and we derive the resulting partial differential equation (pde) 
that the price solves; see equation (2.16). We also present the pde for the price A^'^^ of 
n conditionally independent and identically distributed life insurance risks; see equation 
(2.18). 
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2.1 Mortality Model and Financial Market 

We use the stochastic model of mortality of Milevsky, Promislow, and Young (2005). 
We model the hazard rate for an individual or set of individuals of a given age. If we were 
to model a population's hazard rate, then we would take into account age and time as 
in Lee and Carter (1992) and, more recently, Ballotta and Haberman (2006). However, 
because we consider a single age, we simply model the hazard rate as a stochastic process 
over time. 

We assume that the hazard rate (or force of mortality) of an individual at time 
t follows a diffusion process such that if the process begins at Aq > A for some positive 
constant A, then At > A for all t > 0. From a modeling standpoint, A could represent 
the lowest attainable hazard rate remaining after all causes of death such as accidents 
and homicide have been eliminated; see, for example, Gavrilov and Gavrilova (1991) and 
Olshansky, Carnes, and Cassel (1990). 

Specifically, we assume that 



in which {VF/^} is a standard Brownian motion on a probability space (0,JF, P). The 
volatility a is either identically zero, or it is a continuous function of time t bounded below 
by a positive constant k on [0,T]. The drift /U is a Holder continuous function of A and t 
for which there exists e > such that if < A — A < e, then fi{X, t) > for all t G [0, T]. 
After Lemma 3.2 below, we add more requirements for fi. Note that if cr = 0, then At is 
deterministic, and in this case, we write X{t) to denote the deterministic hazard rate at 
time t. 

Suppose an insurer issues a term life insurance policy to an individual that pays 
$1 at the moment of death if the individual dies before some time T. In Section 2.2, to 
determine the value of the insurance policy, we create a portfolio composed of the obligation 
to underwrite the life insurance, of default-free, zero-coupon bonds that pay $1 at time T, 
regardless of the state of the individual, and money invested in a money market account 
earning at the short rate. Therefore, we require a model for bond prices, and we use a 
model based on the short rate and the bond market's price of risk. 

The dynamics of the short rate {rt}, which is the rate at which the money market 
increases, are given by 



dXt = fi{Xt, t){Xt - X)dt + a{t){Xt - X)dWt^, 



(2.1) 



drt = b{rt,t)dt + d{rt,t)dWt, 



(2.2) 
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in which b and d > are deterministic functions of the short rate and time, and {Wf} is 
a standard Brownian motion with respect to the probabihty space (Q, T ^ P), independent 
of {W^^}. We assume that h and d > are such that > almost surely for all t > 0. 
Assuming that the dynamics of the short rate and the hazard rate are independent is rea- 
sonable because one does not necessarily expect changes in interest rates to affect changes 
in mortality rates and vice versa. However, if one believes that it is more accurate to allow 
the two processes to be correlated, then the work in this paper extends appropriately, 
albeit with added complexity because of the correlaton. 

From the principle of no-arbitrage in the bond market, there is a market price of risk 
process {qt) for the bond that is adapted to the filtration generated by {W^}; see, for 
example, Lamberton and Lapeyre (1996) or Bjork (2004). Moreover, the bond market's 
price of risk at time t is a deterministic function of the short rate and of time, that is, 
qt = q{rt,t). Thus, the time-t price F of a default-free, zero-coupon bond that pays $1 at 
time T is given by 



F(r, t; T) = 



e-I^'^'^'ds 


rt = r 







(2.3) 



in which Q is the probability measure with Radon-Nikodym derivative with respect to P 
given by ^\-p^ = exp Jq q{rs, s)dWs — \ q'^{rs, s)ds^ . It follows that {VF^*^}, with 

= Wt + /q qivs, s)ds, is a standard Brownian motion with respect to Q. 

From Bjork (2004), we know that the T-bond price F solves the following pde: 

Ft + b'^ir, t)Fr + ^d\r, t)Frr -rF = 0, 

F(r,T;T) = l, 

in which = b — qd. Henceforth, we drop the notational dependence of F on T because 
T is fixed (and understood) throughout the paper. 

We can use the pde (2.4) to obtain the dynamics of the T-bond price F{rs,s), in 
which we think of = r as given and s e [t,T]. Indeed, 



dF{rs, s) = {ts F{rs, s) + q{rs, s) d{rs, s) Fr{rs, s))ds + d{rs, s) Fr{rs, s) dWs, 
F{rt,t) = F{r,t). ^^'^^ 

We use (2.5) in the next section when we develop the dynamics of a portfolio containing 
the obligation to pay the life insurance and a certain number of T-bonds. 

2.2 Pricing via the Instantaneous Sharpe Ratio 
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In this section, we first describe our method for pricing in this (incomplete) life in- 
surance market. Then, we fully develop the pricing rule for a single life insurance policy; 
that is, n = 1. Next, we consider some special cases when n = 1. Finally, we present the 
pricing rule for the case n > 1. 

2.2.1 Recipe for valuation 

The market for insurance is incomplete; therefore, there is no unique pricing mech- 
anism. To value contracts in this market, one must assume something about how risk is 
"priced." For example, one could use the principle of equivalent utility (see Zariphopoulou 
(2001) for a review) or consider the set of equivalent martingale measures (Dahl and M0ller, 
2006) to price the risk. We use the instantaneous Sharpe ratio because of its analogy with 
the bond market's price of risk and because of the desirable properties of the resulting 
price; see Section 3. Also, in the limit as the number of contracts approaches infinity, the 
resulting price can be represented as an expectation with respect to an equivalent martin- 
gale measure; see Corollary 4.6 below. Because of these properties, we anticipate that our 
pricing methodology will prove useful in pricing risks in other incomplete markets. 

Our method for pricing in an incomplete market is as follows: 

A. First, define a portfolio composed of two parts: (1) the obligation to underwrite the 
contingent claim (life insurance, in this case), and (2) a self- financing sub-portfolio of 
T-bonds and money market funds to (partially) hedge the contingent claim. 

B. Find the investment strategy in the T-bonds so that the local variance of the total 
portfolio is a minimum. If the market were complete, then one could find an investment 
strategy so that the local variance is zero. However, in an incomplete market, there 
will be residual risk as measured by the local variance. This control of the local 
variance is related to local risk-minimization as described in Schweizer (2001). 

C. Determine the price of the contingent claim so that the instantaneous Sharpe ratio of 
the total portfolio equals a pre-specified value. This amounts to setting the drift of 
the portfolio equal to the short rate times the portfolio value plus the Sharpe ratio 
time the local standard deviation of the portfolio. Bjork and Slinko (2006) apply the 
idea of limiting the instantaneous Sharpe ratio to restrict the possible prices of claims 
in an incomplete market. 

2.2.2 Pricing a single life insurance policy 

Denote the time-t value (price) of a term life insurance policy that pays $1 at the time 
of death of the individual between times t and T by ^ = ^(r. A, t), in which we explicitly 
recognize that the price of the insurance depends on the short rate r and the hazard rate 



6 



A at time t. (As an aside, by writing A to represent the value of the hfe insurance pohcy, 
we mean that the buyer of the pohcy is stiU ahve. If the individual dies before time T, 
then the value of the policy jumps to $1.) 

Suppose the insurer creates a portfolio with value lit at time t. The portfolio contains 
(1) the obligation to underwrite the life insurance of $1, with value —A, and (2) a self- 
financing sub-portfolio of T-bonds and money market funds with value Vt at time t to 
(partially) hedge the risk of the insurance policy. Thus, = —A{rt, Xf,t) + Vt- Let 
TTt denote the number of T-bonds in the self-financing sub-portfolio at time t, with the 
remainder, namely Vt — 7rtF{rt,t), invested in a money market account earning the short 
rate r^. 

The insurance risk cannot be fully hedged because of the randomness inherent in 
the individual living or dying. If the individual dies at time t < T, then the value of 
the insurance policy jumps to $1 because the insurer is obligated to pay $1 immediately; 
therefore, the value of the portfolio 11^ jumps from —A{r, X,t) + Vt to —1 + Vt. In other 
words, the value of the portfolio changes instantly changes by A{r, A, t) — 1, or equivalently, 
the portfolio value decreases by 1 — A{r, A, t), the so-called net amount at risk (Bowers et 
al., 1997). 

We next describe the dynamics of the total portfolio lit by specifying the dynamics 
of its two pieces A{rt, Xt, t) and Vt. By Ito's Lemma (see, for example, Protter (1995), the 
dynamics of the value of the life insurance policy are given as follows: 

1 1 

dA{rt, Xt, t) ^Atdt + Ar dn + -Arr d[r, r]t + Ax dXt + -Axx d[X, X]t + (A - 1) dNt, 

= Atdt + Ar{bdt + ddWt) + ^Arr d^ dt + Ax{Xt - X){iidt + adWt^) 
+ ]^Axx o''{Xt - Xfdt + {A-l) dNt, 

(2.6) 

in which we suppress the dependence of the functions o, 6, etc. on the variable r. A, and t. 
Also, [r, r], for example, denotes the quadratic variation of r, and {Nt} denotes a counting 
process with stochastic parameter At at time t that indicates when the individual dies. 
Because Vt = 7rtF{rt, t) + (Vt — 7rtF{rt, t)) is self-financing, its dynamics are given by 

dVt = TTt dF{rt, t) + n{Vt - TTtFin, t))dt ^{TTtqdFr + n V) dt + 7TtdFr dWt, (2.7) 

in which the second equality follows from (2.5) and we again suppress the dependence of 
the functions on the underlying variables. 
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It follows from equatfons (2.6) and (2.7) that the value of the portfoho at time t + h 
with h > 0, namely Ht+hj equals 



/t+h pt+h 
dA{rsAs.s) + J dVs 



/t+h j-t+h 
V^A{rs, Xs, s)ds + J d{rs, s)(7r,F^(r„ s) - Arir^, A„ s)) dW^ 

ot+h ft+h 



/z-\-n fi-tn 
ais){X, - X)Ax{rs, A„ s) dW^ + (A(r„ A„ s) - l)(diV, - A, ds) 

i-t+h rt+h 

+ / TTsqirs,s)d{rs,s)Fr{rs,s)ds+ / rsUgds, 
Jt Jt 



(2.8) 

in which T>^ is an operator defined on the set of appropriately differentiable functions on 

R+ X (A, oo) X [0, T] by 

V^v ^vt + bvr + \d^ Vrr + l^{X - X)vx + ^cr'^{X - Xf vxx " X{v - 1) - rv. (2.9) 

Note that we adjust dN^ in (2.8) so that Nt — Jq Xg ds is a martingale. 

In this single-life case, the process {11^} is "killed" when the individual dies. If we were 
to consider the price A^'^^ of n conditionally independent and identically distributed lives 
(conditionally independent given the hazard rate), then {Nt} would be a counting process 
with stochastic parameter nXt at time t such that 11^ decreases by 1 — {^A^'^^ — A^'^~^^^ 
when one of the n individuals dies. We consider A^'^^ later in this section and continue 
with the single-life case now. 

We choose tt^ to minimize the local variance of this portfolio. To this end, we calculate 
the expectation and variance of Ht+h conditional on the information available at time t, 
namely J^f First, given lit = H, define a stochastic process Yh for /j, > by 

pt+h i-t+h i-t+h 

y/i = n - / T>^A{rs,Xs,s)ds-\- / TTsq{rs,s) d{rs,s) Fr{rs, s) ds + I Vs^sds. 
Jt Jt Jt 

(2.10) 

Thus, E(nt_|_/i|jFt) = E^''^'*(Y^), in which E^'''^'* denotes the conditional expectation given 
rt — r and At = A. Also, from (2.8) and (2.10), we have 

rt+h 

Tlt+h = Yh+ I d{rs,s){'KsFr{rs,s) - Ar{rs,Xs,s))dWs 

/t+h r-t+h y^-"-^) 

a{s){Xs - X)Ax{rs, A„ s) dW^ + j (A(r„ A„ s) - l)(diV, - A, ds). 
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It follows that 



rt+h 

/t+h pt+h 
a\s){\s - XfAlirs, A„ s)ds + E^'^'* / A,(^(r„ A„ s) - ifds. 

(2.12) 

Recall that we wish to choose tt^ to minimize the local variance lim/j^o Var(nt+/i|.Ft), 
a dynamic measure of risk of the portfolio; therefore, tt^ = ^^.(rt, A^, t)/F^(rt, t). Under 
this assignment, the drift and local variance become, respectively, 

lim i(E(n,+;,|jr,) - n) = -V^" A{r, A, t) + rn, (2.13) 

ft— >o n 

and 



hm -Var(n,+ft|.^,) = a\t){X - Xf Al{r,\t) + KA{r,\t) - l)\ (2.14) 

in which 'D^'^ is given by the expression in (2.9) with h replaced by 6*^ = 6 — qd. 

Now, we come to pricing via the instantaneous Sharpe ratio. Because the minimum 
local variance in (2.14) is positive, the insurer is unable to completely hedge the mortality 
risk underlying the life insurance. Therefore, the price should reimburse the insurer for 
its risk, say, by a constant multiple a of the local standard deviation of the portfolio, a 
measure of the mortality risk inherent in the insurance. It is this a that is the instantaneous 
Sharpe ratio. We could allow a to be a function of say A and t, to parallel the market 
price of risk process {qt} in the bond market. However, we choose a as a constant for 
simplicity. Also, see Milevsky, Promislow, and Young (2007) for further discussion of the 
instantaneous Sharpe ratio. 

To determine the value (price) A of the life insurance, we set the drift of the portfolio 
equal to the short rate times the portfolio plus a times the local standard deviation of the 
portfolio. Thus, from (2.13) and (2.14), we have that A solves the equation 

-V^'^A + rn = rn + a^a^{t){X - XfAl + X{A - 1)^, (2.15) 

for some < a < -\/A- (I^i the proof of Theorem 3.4 the reason for the upper bound on a 
will become apparent.) It follows that A = A{r, A, t) solves the non-linear pde given by 
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At + h'^ir, t)Ar + -d\r, t)Arr + //(A, t){\ - X)Ax + -a\t){\ - XfA^x - rA 

' - XiA - 1) = -a^a''{t){X - X)^Al + X{A - 1)2, (2.16) 
. A(r, A,T) =0. 

Note that if the individual is still alive at time T, then the policy is worthless. Hence, we 
have the terminal condition A{r, A, T) = 0. 

If we had been able to choose the investment strategy {nt} so that the local variance 
in (2.14) were identically zero (that is, if the risk were hedgeable), then the right-hand side 
of the pde in (2.16) would be zero, and we would have a linear differential equation of the 
Black-Scholes type. 

2.2.3 Special cases for n = 1 

If there were no risk loading, that is, if a = 0, then the price is such that the expected 
return on the price at time t is rt + Xt. The rate rt arises from the riskless money market, 
and At arises from the expected release of reserves as individuals die. If a > 0, then 
the expected return on the price is greater than + A^. Therefore, a, the Sharpe ratio, 
measures the degree to which the insurer's total expected return is in excess of rt + Xt, as 
a proportion of the standard deviation of the return. One can think of the right-hand side 
of (2.15) as adding a margin to the return of the portfolio because the insurance contract 
is not completely hedgeable due to the mortality risk. 

Consider the special case for which o" = 0, that is, the hazard rate is deterministic. 
Suppose A('u), for -u > t, is the solution of dX = /u(A,w)(A — A)^^ with initial value 
X{t) = A; then, under the assumption that A > (which we later demonstrate in Theorem 
3.4), (2.16) becomes a linear pde with solution 



A{r,X,t) = F(r,t;s)e-r(^("^+"^)''"(A(s) + av^) ds, (2.17) 

in which F(r, t; s) is the time-t price of a bond that pays $1 at time s >t. The price A in 
(2.17) is similar to the pricing rules encountered in courses in Life Contingencies because 
one can think of A(s) -|- a^/X(s) > as a hazard rate. Note that A in (2.17) increases with 
respect to a, as occurs more generally. 

Below, we show that A increases with a when the hazard rate is stochastic; see 
Theorem 3.7. For this reason, we refer to the solution A of (2.16) as the risk-adjusted 
price for a life insurance contract, in which a controls the risk adjustment. 
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2.2.4 General case of n > 1 

To end this section, we present the pde for the price A^'^^ of n conditionally inde- 
pendent and identically distributed insurance risks. Specifically, we assume that all the 
individuals are of the same age and are subject to the same hazard rate as given in (2.1); 
however, given that hazard rate, the occurrences of death are independent. As discussed 
in the paragraph preceding equation (2.9), when an individual dies, the portfolio value 11 
decreases by 1 — (A^"'^ — A*^"""-*^)), or increases by A^"^^ — A^'^~^^ — 1. By paralleling the 
derivation of (2.16), one can show that A^"^^ solves the non-linear pde given by 

- ryl(^) - nX f - - l) 

, ^ - - (2.18) 

= -a Ja2(t)(A - A)2 (aI^^^)' + nX (AW - A(^-^) - l)^ 
, ^^(r, A,T) = 0. 

The initial value in this recursion is ^(0) = 0, and the price A defined by (2.16) is A(^\ 

3. Qualitative Properties of A^'^^ 

To demonstrate a variety of properties of the risk-adjusted price we need a 

comparison theorem. We begin by stating a relevant one-sided Lipschitz condition along 
with growth conditions. We require that the function g = g{r, X,t,v,p,q) satisfies the 
following one-sided Lipschitz condition: For v > w, 

g{r, A, t,v,p,q)- g(r, A, t, w,p\q')< Ci (r, A, t) {v-w) + C2 (r, A, t)\p-p'\+ C3 (r, \t)\q-q'\, 

(3.1) 

with growth conditions on ci, C2, and C3 given by 

< ci(r,A,t) < K(l + (lnr)2 + (ln(A- A))2), (3.2a) 

< C2(r, A, t) < Kr{l + \ lnr| + | ln(A - A)|), (3.26) 



and 



< C3(r, A, t) < K{X - A)(l + I lnr| + I ln(A - A)|), (3.2c) 
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for some constant K >0, and for all (r, A, t) e R+ x (A, oo) x [0, T]. 

Throughout this paper, we rely on the following useful comparison principle, which 
we obtain from Walter (1970, Section 28). We omit the proof because it is a standard 
application of Walter's work. 

Theorem 3.1. Let G = R""" x (A, oo) x [0, T], and denote by Q the collection of functions 
on G that are twice- differentiate in their first two variables, namely r and A, and once- 
differentiable in their third, namely t. Define a differential operator C on Q by 

1 1 

Cv = vt + -d^{r,t)vrr + -a'^{t){X - Xfvxx+ g{r,X,t,v,Vr,vx), (3.3) 

in which g satisfies (3.1) and (3.2). Suppose v,w & Q are such that there exists a constant 
K >0 with V < e^{(i^'-)^+('^(^-A))'} and w > -e^Wnrf+{H^-y)f} for large (Inr)^ + 
(ln(A — A))^. Also, suppose that there exists a constant K > such that d{r,t) < Kr for 
allr>0 andO<t<T. 

Then, if (a) Cv > Cw on G, and if (b) v{r, A, T) < w{r, A, T) for all r > and A > A, 
then V <w on G. 

As a lemma for results to follow, we show that the differential operator associated 
with our problem satisfies the hypothesis of Theorem 3.1. 

Lemma 3.2. Define gn, for n > 1, by 

gn{r,X,t,v,p,q) = b^ {r,t)p + n{X,t){X - X)q - rv - nX (v - A^""'^^ - l) 

I ^ ('-^^ 

+ a^a^{t)iX - A) V + nX{v- A(^-^) - l) , 

in which ^("■~^) solves (2.18) with n replaced by n — 1. Then, gn satisfies the one-sided 
Lipschitz condition (3.1) on G. Furthermore, if \b^{r^t)\ < Kr{l-\- |lnr|) and |/i(A, t)| < 
K{1 + I ln(A - A)|), then (3.2) holds. 

Proof. Suppose v > w, then 

5fn(r, A, t, V, p, q) - Qnir, A, t, w, p\ q') 

= b^ (r, t){p - p') -\- ii{X,t){X - X){q - q') - {r -\- nX) {v - w) 

+ a^a^{t){X - Xfq^ + nX{v- a("-i) - l)^ 

- a^a^{t){X - A)2(g')2 + nX {w - a^^-^) - l)^ 
< \b^{r,t)\\p-p'\-\-{\i^{X,t)\-\-aa{t)){X-X)\q-q'\ - (r -\- nX - aVnXj {v - w) 
<\bQ(r,t)\\p-p'\ + M\t)\+aait))iX-X)\q-q'\. 

(3.5) 
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In the first inequafity, we use the fact that iiA>B, then y/C^ + -y/C^ + <A-B. 
For the second inequahty, recall that < a < ^/X- Thus, (3.1) holds with ci(r, A, t) = 0, 
C2{r,X,t) = \b^{r,t)\, and C3{r,X,t) = \ii{X,t)\ + aa{t). Note that C2 satisfies (3.26) if 
\b^{r,t)\ < Kr(l + |lnr|), and C3 satisfies (3.2c) if |//(A,t)| < i^(l + | ln(A - A)|). □ 

Assumption 3.3. Henceforth, we assume that the volatility on the short rate d satisfies 
the growth condition in the hypothesis of Theorem 3.1 and that the drifts and satisfy 
the growth conditions in the hypothesis of Lemma 3.2. For later purposes (for example, 
see Theorem 3.6), we also assume that n\ is Holder continuous and satisfies the growth 
condition |/xa|(A - A) + < K (l + (ln(A - A))^). 

For the remainder of Section 3, we study properties of the risk-adjusted price 

^(n) for 

n life insurance contracts. In Section 3.1, we demonstrate two basic properties of A^'^\ 
namely, < ^("^^ < n and A^^^ > 0. In Section 3.2, we examine how A^^^ changes as 
we change the model parameters related to the mortality risk, specifically, a, /i, and a. 
Finally, in Section 3.3, we show that A^'^^ is subadditive with respect to n. 

3.1 Basic Properties of A^'^^ 

In the first application of Theorem 3.1 and Lemma 3.2, we show that < A^"^) < n 
for n > 0. Because the payoff under life insurance is nonnegative, we expect its price to be 
nonnegative. Also, if the hazard rate were arbitrarily large, then we would expect the n 
individuals to die immediately, with the corresponding price for the insurance equal to n. 
In work in Section 4, we sharpen these bounds considerably. Throughout, we encourage 
the casual reader to read the proofs lightly and to focus on the (and sometimes surprising) 
results of the theorems, as well as the corresponding discussions of those results. 

Theorem 3.4. < A^'^^ < n on G, for n > 0. 

Proof. We proceed by induction to prove these bounds. The inequality holds with equality 
when n = because A^^^ = 0. For n > 1, assume that < ^("■~^) < (n — 1), and show 
that < < n. 

Define the differential operator £ on ^ by (3.3) with g = Qn from (3.4). Because A^'^^ 
solves (2.18), we have = 0. Denote by the function that is identically on G; 

then, 

£0 = (nX + aV^^ (^A^"-^) + l) > = CA^^'l (3.6) 
Also, ^(^)(r, A,T) = 0; thus. Theorem 3.1 and Lemma 3.2 imply that < on G. 
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Next, we prove the upper bound. Denote by n the function that is identically n on 
G; then, 



- (nX - aV^) ((n - 1) - yl^"-^)) < = CA^^'l (3.7) 



Additionally, A(")(r, A, T) = < n; thus. Theorem 3.1 and Lemma 3.2 imply that A^"^^ < n 
on G. □ 

We end this subsection by showing that A^^^ > 0. This result is intuitively pleasing 
because we expect the price of life insurance to increase as the hazard rate increases, that 
is, as individuals are more likely to die. We first present a lemma which we use in the 
proof of Theorem 3.6. 

Lemma 3.5. A^") < A^^-^) + 1 onG, for n>l. 

Proof. We proceed by induction. This inequality is true for n = 1 because A^^^ < 1 = 
+ 1 by Theorem 3.4. For n > 2, assume that yl("--i) < ^('^"2) + and show that 

^(n) < ^(n-1) ^ I 

Define a differential operator £ on ^ by (3.3) with g = Qn from (3.4). Because A^'^^ 
solves (2.18), we have LA^'''^ = 0. Also, 



C (yl(^-i) + 1^ = -r + (n - 1)A - - l) + aa{X - X) 

t^a^X - A)2 (a^^'^^Y + {n- 1)X {A(^-^) - A(^-^) - l) 



— a\ 
< = CA^'^\ 



(3.8) 



In addition, ^(")(r, A,T) = < 1 = A^'^'^^r, X,T) + 1. Thus, Theorem 3.1 and Lemma 
3.2 imply that yl^") > + 1 on G. □ 

Lemma 3.5 is interesting in its own right because it confirms our intuition that the 
price of adding an additional policyholder, namely A^'^^ — A^'^~^\ is less than the cost of 
paying the life insurance benefit immediately, namely 1. 

Theorem 3.6. A^^^ > on G, for n > 0. 

Proof. We proceed by induction. We know that the inequality holds when n = because 
^(0) = 0. For n > 1, assume that ^^""^^ > 0, and show that ^ > 0. We apply 
a modified version of Theorem 3.1 to the special case of comparing A^^^ with the zero 
function 0. First, differentiate ^("^'s equation with respect to A to get an equation for 
/(") = A^^\ 
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ft^ + hQf!f^^ + IdV^"^ + (/^a(A - A) + - nA - r)/('^) + nA/^"-^) 
+ + a2)(A - A)/f > + \a\\ - X)'fS - n {a(-^ - A^^-'^ - l) 



= —a- 



a 



\l<y'^{^ - A)2 (/H)^ + nA (ylW - yl(^-i) - l)' 
nA - - l) (/(^) - /(^-i^) 



(3.9) 



^'(A - A)2 (/(-))^ + nA (yl(") - yl(^-i) - l) 
[/(-)(r,A,T) = 0. 
Define a difi^erential operator £ on ^ by (3.3) with g = Qn given by 

^^(A, t, q) = b^p +{f,^{X-X) + fi-nX- r)v + nA/^""^) + (// + a2)(A - X)q 

^a2(A- A)2^;2 + nA (A(-) - A(--i) - l)' 
nA - - 1) - /^""^^) 

I y • 

y (72(A - A)2|;2 + nA - - l)^ 

(3.10) 

Prom Walter (1970, Section 28, pages 213-215), we see that we only need to verify that 
(3.1) and (3.2) hold fov v > = w = p' = q' . 

gn{X, t, V, p, q) - gn{X, t, w, 0, 0) = b^p + {lJ,x{X - X) + fj, - nX - r)v + {fj, + cr'^){X - X)q 
a'^iX - Xy + a^iX - X)\q + nA (yl(^) - A^^-^) - l) ^; 



+ a- 



an 
+ — 



^ct2(A - A)2«;2 + nA (A(-) - A(--i) - if 

(^4(") _ ,4("-i) - 1)^ (yl(") - ^4("-i) _ 1)' 



^ ' V(A - A)2^;2 + nA - - 1)^ JnX (A(^) - - l)^ 



+ ctnA 



^Ja^^iX- A)2^;2 + nA (vK") - vK^-i) - l)^ ^nA (ylW _ ^(n-i) _ if ^ 
< (Ia^aI (A - A) + \fx\ + aa) (^; - 0) + |6^| p + {\fj.\ + + aa){X - X)\q\. 

(3.11) 

Thus, by Assumption 3.3, gn satisfies (3.1) with ci = |//a|(A — A) + +Q!cr, C2 = and 
C3 = (|//| + cr^ + aa){X — A), in which ci, C2, and C3 satisfy the growth conditions in (3.2). 
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Next, note that because = ^ satisfies (3.9), = 0. Also, we have 

£0 = (nX + aVnXj /^""^^ + (n + a/2^/^[/X^ (^("-i) + 1 - ^(")), which is greater than 
or equal to by the induction assumption and by Lemma 3.5. These observations, together 
with /(^) (r. A, T) = 0, imply that = A^^'^ > on G. □ 

3.2 Comparative Statics for A^'^'> 

In this section, we show that as we vary the model parameters, the price A^"^^ responds 
consistently with what we expect. 

Theorem 3.7. Suppose ^ < o-i < 0.2 < a/2; cind let ^("■)''^* he the solution of (2.18) with 
a = ai, for i = 1,2 and n > 0. Then, < on G. 

Proof. We proceed by induction. It is clear that the inequality holds for n = because 
^(o),«i = for i = 1,2. For n > 1, assume that A^''-^)'"! < ^('^-i).«2^ and show that 

^(n),ai ^ ^(n),a2_ 

Define a differential operator £ on ^ by (3.3) with g — Qn from (3.4) with a — ai. 
Because A^'^)'"! solves (2.18) with a = cii, we have = 0. Also, 



|^(72(A-A)2 (|4''^''^^)' + nA (^(n),«2 _^(n-l),ai _1)2 

-y^(T2(A-A)2 (^A^^^'^^^' + nA (^(n),«2 _^(n-l),a2 _1)2| (3.12) 



+ (ai - a2)Y (t2(A - A)2 (^aJ^^'"^) + nA (AH'«2 _ A(--i).«2 _ if 

< (nX - aiV^) (^A^n-l),ai _ ^(n-l),a2^ < q ^ 

In addition, both and satisfy the terminal condition ^('^)'"*(A,T) = 0. 

Thus, Theorem 3.1 and Lemma 3.2 imply that < ^(").«2 on G. □ 

Theorem 3.7 states that as the parameter a increases, the risk- adjusted price 
increases. This result justifies the use of the phrase risk parameter when referring to a. 
We have the following corollary of Theorem 3.7. 

Corollary 3.8. Let f,e the solution of (2.18) with a = 0; then, ^(^)'"0 < ^(")'« for 

all < a < \fX, and we can express the lower hound as follows: = riA'^o^ 
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in which A'^'^ is given by 



Xudu 



A. 



ds. (3.13) 



Proof. It is straightforward to show that nA°'^ solves (2.18) with a = 0, and the result 
follows. □ 

Note that A°'^ is what actuaries call the net premium under the risk-neutral measure 
for stochastic interest and the physical measure for stochastic mortality. Therefore, the 
lower bound of (as a approaches zero) is the same as the lower bound of = 

^(i),q:^ namely, We call ^^^"^ — A"'^ the risk charge per person. At the end of 

Section 4, we show how to decompose this risk charge into the charges for finite portfolio 
and stochastic mortality risks; see equation (4.18). 

Next, we examine how the risk-adjusted price A^'^'i varies with the drift and volatility 
of the stochastic hazard rate. We state the following two theorems without proof because 
their proofs are similar to earlier ones. 

Theorem 3.9. Suppose < ^2 on G, and let A*^""^'^^ denote the solution of (2.18) with 
^ = fji-^ for i^l,2 andn>0. Then, < ^^'^2 on R+ x G. 

From Theorem 3.9, we learn that as the drift of the hazard rate increases, then the 
price of life insurance increases. This occurs for essentially the same reason that the price 
increases with the hazard rate; see Theorem 3.6. 

Theorem 3.10. Suppose < cri(t) < a2{t) on [0,T], and let A^"^'"'* denote the solution 
of (2.18) with a = di, for i = l,2 andn>0. If A^^^'^^ > for all n > 0, or if A^^^'"^ > 
for all n > 0, then < ^n R+ x G. 

From Theorem 3.10, we see that if A^"^^ is convex with respect to A, then the risk- 
adjusted price increases as the volatility of the stochastic hazard rate increases. 

3.3 Subadditivity of A^""^ 

We next show that A^"^^ is subadditive. Specifically, we show that for m, n nonnegative 
integers, the following inequality holds: 

^(m+n) < ^(m) ^ ^(n)_ (^3J4^ 

Subadditivity is a reasonable property because if it did not hold, then buyers of insurance 
could insure risks separately and thereby save money. First, we show that the price per 
risk decreases with n, then we obtain subadditivity as a corollary. 
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Theorem 3.11. ^^^"^ decreases with respect to n > 1 on G. 

Proof. We proceed by induction to show that A^"') < :^^A^'^~^^ for n > 2. For n = 2, 
define a diflFerential operator "D on ^ by (3.3) with g = g2 from (3.4). Because A^"^^ solves 
(2.18) with n = 2, we have VA^^'^ = 0. Also, 



V 



(2^(^)) = 2a I ^a^{X-Xr (a^) V 2A QaW - 0' 

I 2 1 (3-15) 

-Ja^{X - A)2 (a^^^^ + (1)A {AW - AW - l)^ i 

<0, 

and A(2)(r, A,T) = 2A(i)(r, A, T). Thus, it follows from Theorem 3.1 and Lemma 3.2 
that < 2^(1) on G. 

Now, suppose for some n > 3, A^^-^) < ^A^'"-^), and show that A^""^ < ^A^''-^\ 
Define a differential operator £ on ^ by (3.3) with g = Qn from (3.4). Because A^"^^ solves 
(2.18), we have CA^'^^ = 0. Also, 



n — 1 I \n— 1 



<n(X- a^) - A("-2)^ < 0. 

The first inequality in (3.16) is not obvious, so we prove it now. Let B denote 
(t(A — A)yl^"' then, the inequality is equivalent to 



(3.16) 



Ib^ + hX Q^^""'^ - < VA (^(n - 1)71(^^-2) - (n - 2) A^^^-^)) 

+ ^JB^ + {n- 1)A - Ai^-^) - if. 

Square both sides of (3.17), cancel B'^ from both sides, and note that the resulting inequal- 
ity holds for all B if and only if the following holds 
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r/ X / m2 (3-18) 

Take the square root of both sides and simphfy to obtain the equivalent inequahty 



< ((n - 1) + V^^) + A^""-^^ (^-{n - 2) - V^T^ + -^^ + VrTH: 



n — 1 
(3.19) 

By the induction assumption, we have that < n=i^(n-2)_ rj,^^^ ^^g -^g^ j^^^^^ .j^j^^ 

following inequality holds: 







< !^^(n-i) ((^ _ 1) + V^^^Ti) + ^(n-1) (r^ - 2) - 

, ^ (3.20) 

+ Vn-1 



which simplifies to 

< 



vn— T — 



n — 1 



1 —A^'^-'^A . (3.21) 

n — 1 I 



By Theorem 3.4, we know that (3.21) holds; thus, we have shown the first inequality in 
(3.16). 

Because C [^A^""-^)^ < = and ^(")(r, A,T) = = ^A(^-i)(r, A, T), it 

follows from Theorem 3.1 and Lemma 3.2 that < on G. □ 

— n— 1 

Corollary 3.12. If m and n are nonnegative integers, then ^('^+"') < A^'^^ + A^'^^ on G. 
Proof. By Theorem 3.11, < and < A^"). Add these two 

' m+n — m+n — 

inequalities to obtain inequality (3.14). □ 

In the next section lies our main result. We show that the limit (as n goes to infinity) 
of ■^A^'^^ solves a linear pde, and we give that pde. 

4. Limit Result for Pricing n Life Insurance Contracts 

In this section, we answer the question motivated by Theorem 3.11; that is, we deter- 
mine the limiting value of the decreasing sequence ^A^"^^ and show that the limiting value 
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solves a linear pde, a rather surprising result. First, we show that ^^^"^ is bounded below 
by P, in which P solves 

Pt + 6^(r, t)Pr + \d\r, t)Prr + (//(A, t) + aa{t)){X - X)Px + ^a\t){X - XfPxx 
-rP- A(P- 1) = 0, 

[P{r,X,T) = 0. 

(4.1) 

Intuitively, P is less than ■^A^'^^ because we replaced the square root in ^^'^^'s pde with 
the square root of the first term. Later, we show that ^^'■"^ equals P as n approaches 
infinity. For this reason, we can view a as the market price of mortality risk (in the limit). 

Theorem 4.1. nP < A^^) on G, for n > 0. 

Proof. We proceed by induction. The inequality holds for n = because A^'^^ = 0. For 
n > 1, assume that (n — 1)P < A^'^~^\ and show that nP < A^'^\ Define a differential 
operator £ on ^ by (3.3) with g = Qn given by (3.4). Because A^"^^ solves (2.18), CA^"^^ — 0. 
Also, 



C{nP) = nX (^A(^-i) - (n - 1)P)) 



+ na}^^a\X-X)^Pl+^{nP-A(r^-^) -if -aiX-X)Px^ (4-2) 
> = £A^''\ 

In addition, nP(r, A, T) = A^'^^r, A, T) = 0; thus. Theorem 3.1 and Lemma 3.2 imply that 
nP < yl(^) on G. □ 

Next, we show that A^'^^ < B^^\ in which B^"^^ solves 



5f ) + 6^(r, t)5^") + -(fir, t)Bi'^^ + (//(A, t) + aa{t)){X - X)B^^^ 

+ ^a\t){X - XfB^;;^ - tB^^^ - (nX + aV^) (s^^) - A^^-^^ - l) = 0, ("^-3) 
[5(-)(r,A,T) = 0, 

in which B^'^^ = 0. Intuitively, A^"^^ < B^'^^ because we replaced the square root in 
pde with the sum of the square roots of the two terms. Recall that \/a^ + 6^ < |a| + 

Theorem 4.2. yl^") < B^") on G, forn>0. 
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Proof. We proceed by induction. The inequality holds for n = because both and 
= are identically zero. For n > 1, ^(^-i) < S^^-^), and show that < S^^). 
Define a differential operator £ on ^ by (3.3) with g = hn, in which hn is given by 

hn{r, A, t, v,p, q) = b^{r, t)p+ (/u(A,t) + aa{t)){X - X)q - rv 

- (nX + aV^) (v - - l) . ^^'^^ 

It is straightforward to show that under Assumption 3.3, ^ = /in in (4.4) satisfies (3.1) and 
(3.2). Therefore, we can apply Theorem 3.1. Because 5^") solves (4.3), £5^^) = 0. Also, 

= aa(X - X)A^;^^ + ay/X (yl^^-^^ + 1 - yl^^)) 



- ^c72(A - xy (4"^)' + A (^(-) - ^(--1) - 1)' 
>o = >cs(^). 

Because both A^"^^ and S^"^-* are zero when t — T, Theorem 3.1 implies that A^"^^ < B^"^^ 
on G. □ 

Finally, we come to the main result of this paper, namely, that limn^oo 77 

Theorem 4.3. 

lim -yl^'^) = P onG, (4.6) 

n— >oo 77, 

in which P solves (4.1). 

Proof. By Theorems 4.1 and 4.2, the theorem is proved if we show that -^B^"^^ — P goes 
to zero as n goes to infinity because — P > — i^*^") > 0. 

Define ^^"^^ on G by = -^B^^^ — P, so the theorem is proved if we show that 
limn_^<x) ^^^^ =0 on G. The function ^'^'^^ solves the pde 

' + 6Q(r, 0$^") + ^ci2(r, t)^t^ + (/.(A, t) + aa{t)){X - A)*^:) 
+ ^a2(t)(A - A)2$J^^) - - (nX + av^) 

= ((n - 1)A + «v7a) P-[x^ asj^ A^^'^^ - 
^$('^^(A,T) = 0, 

with < = 5(1) - P < 1 because S^*^) < n similarly to Theorem 3.4. 

We next show that if $("~^) < K^-i with K^-i a nonnegative constant, then ^^"^^ < 
Kn via a comparison argument with a differential operator based on the pde in (4.7), in 
which 



(4.7) 
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1 n — 1 

Set = 1 because we know that ^^^^ < 1. To that end, define £ on ^ by (3.3) with g 
given by 



5f(r, A, f g) = 6'^(r, t)p + (//(A, + Q;(j(t))(A — X)q — rv — ^nA + a\fnX^ v 

- (^(n- l)A + ay^) P+ ^A + a^^j +ay^. ^ ^ 

It is straightforward to show that under Assumption 3.3, g in (4.9) satisfies (3.1) and (3.2). 
Therefore, we can apply Theorem 3.1. Because ^'^"^^ solves (4.7), we have = 0. 

Assume that < K^-i and show that ^^"^^ < K^- Denote the function that is 

identically in (4.8) by Kn; then. 



'A 



-rKn - (nX + a^^^ Kn - {{n - 1)A + ayfnX^ P+i^ + ^\[^ ^^""^^ 



n 



= -rKn - (nX + aV^) + (nA + aV^) (^T^'^^''"^^ ~ ^) 

+ c./f(l-P). 

(4.10) 

Because ^^(^-i) - P < - P = < it follows from (4.8) and 

(4.10) that 



-^{nX + aVnX) + a^j^ (1 - P) 
[nX + + a^j^ < = 



(4.11) 

1 



< -rKn - 

in which the second inequality follows from < P and the third from a < < on G. 
We also have $(")(r. A, T) = < Kn\ thus, from Theorem 3.1, we conclude that < Kn 
on G. 

Define L„ = nKn- Recall that = 1; thus. 
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Ln = Ln-i + n>2, (4.12) 

from which it foUows that 

""1 dx 

L„ = l + y— <1+/ —<1 + 2V^, n > 2. (4.13) 

FinaUy, we have 

< i^n < - + 4=' ^ > 1' (4-14) 

The right-hand side of inequality (4.14) goes to zero as n goes to oo; thus, goes to 
zero as n goes to oo. □ 

We have the following corollaries of Theorem 4.3. If the hazard rate is deterministic, 
then in the limit, there is no mortality risk. The first corollary demonstrates that the price 
in our setting observes this desirable property and converges to the net premium. 

Corollary 4.4. If a = 0, then lim„^oo ^^^"^ = ^"'^ on G, in which A°''^ is the net 
premium given in (3.13). 

Proof. From Theorem 4.3, we know that ■^A^'^') goes to P as n goes to infinity. Therefore, 
the corollary follows because P — A"^ when cr = 0, which is clear from (4.1) by setting a 
equal to 0. □ 

On the other hand, if the hazard rate is truly stochastic, then in the limit, mortality 
risk remains. The next corollary shows that the limiting price in this case is strictly greater 
than the net premium when t <T. 

Corollary 4.5. If a is uniformly bounded below by n > 0, then lim„^oo -^^^^^ > on 
G, with equality only when t = T. 

Proof. This result follows from the fact that P > A"^, with equality only when t = T. 
Indeed, define a differential operator £ on Q hy (3.3) with = ^fi in (3.4) and with a = 0. 
Thus, = 0, and 

LP = -aa{X - X)Px < = (4.15) 

Also, ^"0(r, A, T) = P{r, A, T) = 0; thus. Theorem 3.1 and Lemma 3.2 imply that P > 

To show that P < for t e [0, T), consider the pde of C = P - The function 
C solves 
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Ct + h^Cr + -(fCrr + (/^ + "(t) (A - A)Ca + CKCT^a 
1 



aO 



+ ^<7'(A - A)2Caa - rC - a (C - 1) = 0, 
tC(r,A,T) = 0. 

From the linear pde in (4.16) and the Feynman-Kac Theorem, we deduce that 



(4.16) 



C(r, A, t) 



eQ) 



{Xs + aa{s)Al\rs, A„ .)) e" K ^^^+^^)'^- ds 



(4.17) 



A 



in which {Ag} follows the diflFusion dXg = {lJi{Xs, s) + aa{s)){\s — X)ds + a{s){Xs — X)dW. 
with Wg = Wg—as. Also, {r^} follows the diffusion dvg = h^{rs, s)ds-\-d{rs, s)dW^ , with 

= Wg + /q* ?(r^, w)(iw; see the discussion following (2.3). Note that a is analogous to 
the bond market's price of risk {g(rs, s)}. The processes {M^g^} and {Wg'^} are independent 
standard Brownian motions with respect to a suitably-defined probability space, and 
denotes expectation on that space. 

From the representation of C = P - in (4.17), we see that C(r, A, t) > for t < T 
because A^ > A and > 0. □ 

After all this work, we can finally decompose the risk charge ■^A'^'^^ — - first 
mentioned in the Introduction and again following Corollary 3.8 - into its component 
risk charges: one for a finite portfolio and another for stochastic mortality. Recall from 
Theorem 4.3 that lim„^oo ^-4(") = P, in which P solves (4.1). Therefore, define ^A^""^ -P 
as the risk charge (per risk) for holding a finite portfolio, and define P — A"^^ as the risk 
charge for stochastic mortality even after selling to an arbitrarily large group. Thus, we 
have 



l^(n) _ ^aO ^ ( l^(n) _ p\ ^ (p _ j^aQ\ ^ /4jg) 

n \n J 

in which the risk charge for stochastic mortality, namely P — is zero if u = by 
Corollary 4.4 and is positive (for t < T) if cr > k > by Corollary 4.5. 

We end this section by more closely examining the limiting value of ^A^'^\ namely 
P, and show that it can be represented as an expectation with respect to an equivalent 
martingale measure. Note the similarity between the limiting price P in (4.1) and the net 
premium in (3.13), or equivalently (2.16) with a = 0. 
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Corollary 4.6. One can represent the limiting value of -A^'^\ namely P, as follows: 



'(r, A, t) = F{r, t; s) E^'* A.e" T ^"'^^ 



ds, 



(4.19) 



in which E is the expectation on the probability space {Q,J^,P), in which the Radon- 
Nikodym derivative of P with respect to P is given by 



dP 
dP 



= exp ( aWt ~ -a t ] . 



(4.20) 



Note that {W^} as defined in the paragraph following (4.17) is a standard Brownian motion 
on {fl.,J-',P, and recall that {Xs} has dynamics 



dXs = {i^{Xs, s) + aa{s)){Xs - X)ds + a{s){Xs - X)dW^. 



(4.21) 



Proof. From the Feynman-Kac Theorem apphed to P's hnear pde (4.1); see Karatzas and 
Shreve (1991), we can represent P as foUows 



P(r,A,t)= (eQ)''^'* j\,e- r^'-^^-'^''^ 



ds 



(eQ) 



-i 



F{r, t; s) E 



x,t 



XsC 



.-I'X^du 



ds 



ds. 



(4.22) 



Note that the representation of P in (4.19) is quite similar to a net single premium 
as defined in courses on Life Contingencies. Indeed, F is the monetary discount function. 



and/(A,t;s) = E^'* 



XgC 



I Xudu 
•J t 



plays the role of the probability density function of the 

time of death. Also, note that we can see clearly in (4.21) that a in the insurance market 
is parallel to the bond market's price of risk {q{rs, s)}. 

It is clear from Corollary 4.6, for example, that we can essentially use any (reasonable) 
model for bond pricing. Additionally, one could use a different sort of default-free bond, 
such as a consol bond (or perpetuity). The latter would be more appropriate in the setting 
of whole life insurance. Alternatively, an insurer could continually roll money into longer- 
term bonds as the shorter-term bonds mature. 

For later reference, note that / solves the pde 
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(ft+ ifiiX, t) + aait))iX - A)/a + ^a^t^X - Xfhx - A/ = 0, 
\/(A,s;s) = A. 
In the next section, we examine a numerical example. 



(4.23) 



5. Numerical Example 

We want to determine how well the lower bound P in (4.1) approximates the price 
A in (2.16), or more generally ^A^^^ in (2.18), for a specific numerical example. Because 
P is a lower bound, then we know that it is biased low. PYom (4.14), it follows that 
P < < iS^'^) < P + i + Due to the slow rate of convergence of -S^^) (and 

presumably of ^A^'^'>) to P (namely, of the order 0{1/ y/n)), we also consider the upper 
bound B — B^^^ that solves (4.3) with n = 1. If S provides a better estimate for A than 
P, then it will provide a better estimate for for n small enough. 

The upper bound B solves 

Bt + b^{r, t)Br + ^(f{r, t)Brr + (/^(A, t) + aa{t)){X - X)Bx 

< +lcT\t){X-XfBxx-rB-(nX + aV^^ (5 - 1) = 0, ^^''^^ 
^i?(r,A,T) = 0, 
Similarly to Corollary 4.6, we can represent B as follows: 

B{r,X,t) = F{r,t;s)E^'' 

Let g denote the function given by ^(A, t; 
g solves 




(gt + (/i(A, t) + aait))iX - X)gx + lcT\t)iX - Xfgxx - (A + a^)g = 0, 
[ g{X, s;s) — X + a\/X. 

Note the parallel between g in (5.3) and / in (4.23). 

It follows from (4.19) and (5.2) that to compute P or B numerically, it is enough to 
compute / or g, respectively, and then numerically integrate the result with F. Because 
the non-linearity in the equation for A arises from the stochastic hazard rate, to test the 
tightness of the bounds P and S, we assume that the interest rate r is identically zero 
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and, thereby, focus on the effect of the stochastic hazard rate. Note that in this case, the 
discount function F is identically 1. 

For this example, we assume that the hazard rate follows the diffusion in (2.1) with 
jj, and a constant. In this case. 

At = A + (Ao - A) exp - ^a^^ + aW^^^ , (5.4) 

a stochastic version of Makeham's law (Bowers et al., 1997). 

See Appendix A for an algorithm for numerically computing the lower bound P. The 
one for the upper bound B is obtained similarly, and for the sake of space, we omit it. 
Because we assume that r is identically zero, we can also solve for the risk-adjusted prive 
A numerically by using a similar algorithm, and this allows us to see which of P or i? gives 
us the better approximation to A. See Appendix B for an algorithm to compute A. 

We use the following values of the parameters for the computation in Table 1: 

• The term of the insurance is T = 10 years. 

• The minimum value of the hazard rate is A = 0.02. 

• The drift of the hazard rate fj, = 0.04. 

• The volatility of the hazard rate cr = 0.10. 

• The risk parameter a = 0.10. 

In Table 1, for a variety of initial values of the hazard rate, we present the limiting 
price P, the risk-adjusted price A, and the upper bound B. For reference, we also include 
the physical probability of dying A""^. Note that in each case, the upper bound B provides 
a better approximation to the risk-adjusted price A than does the limiting price P. For 
other numerical experiments (not shown here), we found the same phenomenon, namely 
that P is a better approximation to A than P. 



Ao 




P 


,4 


D 


0.020 


0.1813 


0.1817 


0.2896 


0.2897 


0.021 


0.1914 


0.1919 


0.3010 


0.3017 


0.022 


0.2014 


0.2025 


0.3126 


0.3139 


0.023 


0.2112 


0.2128 


0.3237 


0.3256 


0.024 


0.2214 


0.2235 


0.3352 


0.3377 


0.025 


0.2300 


0.2326 


0.3449 


0.3477 


0.030 


0.2763 


0.2812 


0.3953 


0.4004 


0.035 


0.3187 


0.3256 


0.4397 


0.4466 


0.040 


0.3609 


0.3696 


0.4826 


0.4909 


0.050 


0.4338 


0.4451 


0.5536 


0.5639 
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0.060 0.5017 0.5150 0.6169 0.6285 

0.070 0.5530 0.5675 0.6630 0.6753 

6. Summary and Conclusions 

We developed a risk-adjusted pricing method for life insurance by assuming that 
the insurance company is compensated for its risk in the form of a given instantaneous 
Sharpe ratio of an appropriately-defined (partially) hedging portfolio. Because the market 
for insurance is incomplete, one cannot assert that there is a unique price. However, 
we believe that the price that our method produces is a valid one because of the many 
desirable properties that it satisfies. In particular, we studied properties of the price of n 
conditionally independent and identically distributed life insurance contracts. In Theorem 

3.11, we showed that the risk charge per person decreases as n increases, and in Corollary 

3.12, we showed that the price is subadditive with respect to n. 

Arguably our main results are those dealing with the limiting price per person. Our 
most important result is that the limiting price solves a linear pde (Theorem 4.3), and we 
provided a probabilistic interpretation of the limiting price as an expectation with respect 
to an equivalent martingale measure (Corollary 4.6). Our work, therefore, extends that of 
Blanchet-Scalliet, El Karoui, and Martellini (2005) and Dahl and M0ller (2006). 

We proved that if the hazard rate is deterministic, then the risk charge per person goes 
to zero as n goes to infinity (Theorem 4.3 and Corollary 4.4). Moreover, we proved that if 
the hazard rate is stochastic, then the risk charge person is positive as n goes to infinity, 
which reflects the fact that the mortality risk is not diversifiable in this case (Theorem 
4.3 and Corollary 4.5). Additionally, in equation (4.18), we decomposed the per-risk risk 
charge into the finite portfolio and stochastic mortality risk charges. 

Milevsky, Promislow, and Young (2005, 2007) proved similar properties for the risk- 
adjusted prices of pure endowments and annuities; see Milevsky, Promislow, and Young 
(2007) for an elementary discussion of this general topic. Because of these properties, 
we anticipate that our pricing methodology will prove useful in pricing risks for many 
insurance products. For example, Bayraktar and Young (2007a) apply the method in this 
paper to price both pure endowments and life insurance. They showed that the price of the 
two products combined is less than the sum of the individual prices, an intuitively pleasing 
result. An interesting extension would be to allow the hazard rate to exhibit jumps, as a 
model for catastrophes. 

We also believe that our valuation method will be useful in pricing risks in other 
incomplete markets. For example, Bayraktar and Young (2007b) price options on non- 
traded assets for which there is a traded asset that is correlated to the non-traded asset. 
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Then, they to price options in the presence of stochastic volatihty, and the instantaneous 
Sharpe ratio in this case equals the market price of volatihty risk, similar to interpreting 
the instantaneous Sharpe ratio in this paper as the market price of mortality risk in (4.1). 



Appendix A. Algorithm for computing P 

We present this algorithm in some detail so that the interested reader can reproduce 
our results (or related results). First, we transform the equation for /, namely (4.23), by 
defining t = s — t, y — ln(A — A), and /(?/, r) = /(A + e^, s — r). Then, / solves 

/r=A/. + ^^'/..-(A + e^)/, ^^^^ 
/(y,0) = A + e^ 

in which (1 = ji + aa — \(t'^. Instead of solving (A.l) for (j/, r) in the domain R x [0, T] , we 
solve (A.l) on the domain [— M , M ] x [0,T], for M such that is approximately zero. 
Therefore, we require boundary conditions f{—M,T) and f{M,T). 

If X^ = A, then Ag = A for all s > t. Thus, /(A, t; s) = Ae~^(*~*), and it follows that 
the appropriate boundary condition at y = —M is f{—M,T) = Xe~-'^. 

If Xt is large, then we expect the individual to die immediately, so that /(A, s) is 
approximately for s > t. Thus, the boundary condition at y = M is f{M,T) = for 
T > 0. 

Subdivide [— M , M] x [0, T] into a grid of rectangles of size h by k. Thus, define 
= -M + n/i forn = 0, 1, . . . , TV with N = 2M/h, and define tj = jk for j = 0, 1, . . . , J 
with J = T/k. Assume that M and h are such that N is an integer; similarly, assume that 
k is such that J is an integer. Now, define fnj = fiUm'^j)- 

We use a backward difference (or implicit) scheme because it is inherently more stable 
(DuChateau and Zachmann, 1986). Note that 

friyn.Tj) = ^^^±P^ +0(/c), (A.2) 

fviVn, rj) = + 0{h'), (A.3) 



and 



f fn, ^\ /n+l,j+l '^fn,j+l + fn-l,j + l r\( u'i\ (\ A\ 

hyiVn, Tj) = + 0[h ), (A.4) 



so that the following equation approximates (A.l) at (y, r) = {yn, Tj) up to order 0{k+h'^) 
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fn,j + l — fnj _ ^ fn+l,j + l — fn-l,j + l 1 2 /^+1,3 + 1 ~ ^fn,j + l + fn-l,j + l 

k ~^ 2h ^2 (A.5) 



If we define a = p,^ - ^cr^, 6=1 +cr^ + kX, and c = ~ 7^' ^^^^ (^-5) 

becomes 



afn-i,j+i + {b + key")fnj+i + cf^+ij+i = fnj-, (A.6) 
for n = 2, 3, . . . , A?" — 2 and j = 0, 1, . . . , J — 1, with the boundary conditions giving us 



(6 + ) h,3+i + c/2,,+1 = hj - aAe-^(^■+l)^ for j = 0, 1, . . . , J - 1, 



(A.7) 



and 



a/AT-ij+i + (6 + A;e2'--) /at-u+i = /iv-i,„ for j = 0, 1, . . . , J - 1. (A.8) 
The system of equations in (A.6) - (A.8) can be represented in matrix form by 



Mf,+i = f,- 



aAe-^^^+^^^ 0, . . . , 0] * , for j = 0, 1, . . . , J - 1. 



(A.9) 



in which the superscript t denotes matrix transform and M is the tri-diagonal matrix with 
main diagonal given by the terms h + ke^^ , h + ke^^ , . . . , 6 + ke^^-^ , with the sub-diagonal 
identically a, and with the super-diagonal identically c. Finally, fj denotes the column 

r _ _ _ -1 t 

vector /ij, /2j, ...Jn-i,j , and similarly for f^+i. 

Beginning with the initial value f^o = A + e^" for n = 1, 2, . . . , — 1, we solve (A.9) 
repeatedly until we obtain /„j for ?i = l,2,...,A — 1. We use the Thomas algorithm for 
solving (A.9), which relies on the LU-decomposition of M into a product of a lower- and 
an upper-diagonal matrix. 

Once we have /, then we want to estimate P at time 0. 



'(r,A,0)= / F{r,0;s)f{X,0;s)d. 
Jo 

^F(r, 0; 0)/(A, 0; 0) + F{r, 0; jk)f{X, 0; jk) + ^^(r, 0; T)/(A, 0; T) 



j-i 



k. 
(A.IO) 
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Suppose that A = A + e^ = A + e^", then (A. 10) becomes 



P(r,A,0) 



j-i ^ 



(A.ll) 



which becomes the following simplified expression because we assume that r is identically 
zero, or equivalently F is identically one: 



P(r,A,0) 



1 



■^-^ ^ 1 _ 

/nO ~l" ^ ^ /nj ~l" '^fnJ 



k. 



(A.12) 



In a sense, the expression in (A.12) is a modified probability of dying before time T. It will 
be larger than the physical probability of dying, which is measured when a = 0. However, 
P will be less than the corresponding value for the price A and the latter less than the 
upper bound B. 

In pseudo-code form, here is the detailed algorithm for computing f^j for n = 1,2,..., 
N-1 and for j = 1,2, . . . , J: 

1. Define M so that is approximately zero, or equivalently, so that is quite large. 

2. Define T as the horizon for the life insurance policy. To price whole life insurance, set 
T equal to some large value, such as 300. 

3. Define the ^/-step-size h, so that N = 2M/h is an integer. 

4. Define the r-step-size k, so that J = T/k is an integer. 

5. Define the parameters of the stochastic hazard diffusion n, a, and A, and define the 
pricing parameter a e [0, V^)- Then, define fi = n + aa — \(t'^. 

6. Define a = fi^ - ^cj^^, b=l + a^^ + kX, and c = -fi^ - i-^ A 



2ft 2*^ ft2' 



7. Forn = 1, 2, . . . , iV - 1, define yn = -M + nh. 

8. Forn = 1, 2, . . . , iV - 1, define fno = \ + e"^ . 

9. Define h = h + ke^^ . Then, forn = 2, 3, . . . , iV - 1, define In = b + ke^^ - ac//,,_i. 
10. For j = 0, 1, . . . , J- 1, define zij = i (^fij - aAe-^(j+^)'=) , and for n = 2, 3, . . . , iV- 1, 

define Znj = j- ifnj — «^n-i,j)- Then, still within the j-loop, define /tv-ij+i = 



ZN-i,j and for n = TV - 2, TV - 3, . . . , 1, define fn,j+i = Znj - ^fn+i,j+i- 



For the work in this paper, we set M 
22,026), h = 0.1, and k = 0.01. 



10 (e 



-10 



4.54 X 10 ^, or equivalently, e 
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Appendix B. Algorithm for computing A 
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The algorithm to compute A is based on the one in Appendix A, so we simply write 
the algorithm. One difference between this algorithm and the one in Appendix A is that 
we approximate the derivative Ax in the non- linear term on the right-hand side of (2.16) 
with a forward difference and include that non-linear term on the right-hand side of the 
analog of (A. 9). Another difference lies in the boundary conditions for small and large 
values of A: A{0,X, t) = 1 - e-U+«V^)(^-*), and A{0, X,t) = l for A large. 

In pseudo-code form, here is the detailed algorithm for computing A^j forn = 1,2, ... , 
N — 1 and for j = 1, 2, . . . , J when r is identically zero: 

1. Define M so that is approximately zero, or equivalently, so that is quite large. 

2. Define T as the horizon for the life insurance policy. To price whole life insurance, set 
T equal to some large value, such as 300. 

3. Define the j/-step-size h, so that N = 2M/h is an integer. 

4. Define the r-step-size k, so that J = T/k is an integer. 

5. Define the parameters of the stochastic hazard diffusion /i, a, and A, and define the 
pricing parameter a e [0, Then, define jj! = jj. — |cr^. 

6. Define a = ^i'^- \a^^, 6 = 1 + + /cA, c = - h^W^ ^nd G = ^a^ {^)\ 

7. Forn = 1, 2, . . . , TV - 1, define yn = -M + nh. 

8. Forn = 1,2, ...,iV- 1, define ino = 0. 

9. Define h = b + ke^^ . Then, for n = 2, 3, . . . , iV - 1, define ln = b + keV^ - ac/ln-i- 
10. For j = 0,1,..., J- 1, define 



1 



V2 2 
G [A2j - 1 + e-U+"VA)^'=) + k2 (A + efi) (^ii, - l) , and for 



n = 2, 3, . . . , iV - 2, define 



Znj — , 



Anj + k{X + e^") - azn-i,j + aS. 



nj 



in which Snj — \j G j — An-i,j^ -|- /c^ (A + e^^") i^A^j — 1 j . Finally, define 



ZN-l,j - -J 

I'M-! 



^N-ij + /c(A-|-e^^-') - c - azN-2,j +aSN-i,j 



I 2 2 

in which Snj = J G (l - AN-2,j^ + k'^ {X + eu^-^) (An-i,] - l) • Then, still with- 
in the j-loop, define Ajv-ij+i = ^N-ij and for n = N — 2, N — 3, . . . ,1, define 
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